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Abstract 

The new BABAR data on the pion-photon transition form factor arouses people's new interests 

^ ; 

♦^1 on the determination of pion distribution amplitude. To explain the data, we take both the leading 

Q ' valence quark state's and the non-valence quark states' contributions into consideration, where the 



valence quark part up to next-to-leading order is presented and the non- valence quark part is esti- 
mated by a phenomenological model based on its limiting behavior at both Q^ — )• and Q^ — )• oo. 



(N 
> 

jy^ I Our results show that to be consistent with the new BABAR data at large Q^ region, a broader 

ff^ ■ other than the asymptotic-like pion distribution amplitude should be adopted. The broadness of 



the pion distribution amplitude is controlled by a parameter B. It has been found that the new 
BABAR data at low and high energy regions can be explained simultaneously by setting B to 
be around 0.60, in which the pion distribution amplitude is closed to the Chernyak-Zhitnitsky form. 
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I. INTRODUCTION 

The pion-photon transition form factor 77* — )■ vr*^, which relates two photons with one 
hghtest meson, is the simplest example for the perturbative application to exclusive pro- 
cesses. It provides a good platform to study the property of pion distribution amplitude 
(DA), i.e. one can extract useful information on the shape of the leading- twist pion DA by 
comparing the estimated result on the transition form factor F^^((5^) with the measured 
one. The CELLO collaboration has measured the pion-photon transition form factor a long 
time ago, where, one of the photons is nearly on-shell and the other one is off-shell with a 
virtuality in the range of low energy region (Q"^ < 3 GeV^) jj. Later on, the CLEO collab- 
oration also measured such form factor but with a broader range of Q^ E [1.5, 9.2] GeV^ j^]- 
Very recently, BABAR collaboration does a more precise measurement at both low and high 
energy region, and their data shows that in the range of Q^ G [4, 40] GeV^, the pion-photon 
transition form factor behaves as [3| 

where A = 0.182 ± 0.002 and B = 0.25 ± 0.02. Such large Q'^ behavior contradicts the well- 
known asymptotic prediction j^, i.e. Q'^Ft,^{Q'^) tends to be a constant (2/7r) for asymptotic 



a. 



DA (pasix, Q )\q2-^oo = 6x(l — x), where the pion decay constant /,r = 92.4 ± 0.25 MeV 
By extending the previous next-to-leading order (NLO) corrections [g, |7| to the present 
large Q^ region, or even by including the next-to- next-to-leading order corrections |8|, l9|] , the 
significant growth of the pion-photon transition form factor between 10 and 40 GeV^ cannot 
be explained by using the asymptotic or asymptotic-like DA. 

Therefore, many attempts have been tried to solve the present puzzle, some authors have 
been argued that in contrary to the conventional adopted asymptotic- like DA, the pion DA 



should be quite broad or even fiat in its whole region 10l-ll3l|. More explicitly, with a flat 
DA (p(x) = 1, Ref. [lO(] shows that the present BABAR data at large Q^ can be explained by 
choosing proper values for the phenomenological parameters for the logarithmic model and 
the Gaussian model constructed there. However, there is no strong reason to support such 
a flat DA, since the introduced infrared regulator m"^ (or a) is rightly fitted by the BABAR 
data. Moreover, one may observe that Ref. [lO|] fails to explain the small Q^-behavior, and 
it can not reproduce the well-known value of Ft^^{Q^ = 0) = 1/(47?^/,^) that is derived from 



measuring the rate of vr*^ — )■ 77 [ij]. Also it can be easily seen that the fiat DA with the 
wavefunction model suggested in Ref. [lO| can not derive the right behavior at Q^ — )■ 0, since 
as will be shown later it will lead to the probability of finding the valence quark state in the 
pion, Pqq = Jq (g^^j^^) dx, and the charged mean squared radius, (r^+)59 = j^ (^^) dx, 
both of which are divergent. Furthermore, with such a fiat DA, the end-point singularity 
shall be emerged in many exclusive processes, such as 5 — > light meson transition form 
factors, which malies them not calculable in perturbative QCD. This shall greatly compress 
the applicability of perturbative QCD ^. 

At present, there is no definite conclusion on whether pion DA is in asymptotic form 
1^, in Chernyak-Zhitnitslcy (CZ) form 17| or even in fiat form [18!]. The pion DA can be 



expressed in Gegenbauer expansion ^ The value of the Gegenbauer moments have been 
studied in various processes, cf . Refs. 19l-l27j . The lattice result of Ref. |27| prefers a narrower 



DA with 02(1 GeV^) = 0.07(1), while the lattice results 25|, l26| prefer broader DA, i.e. they 



obtain 02(1 GeV^) = 0.38 ± 0.23T^:o6 and 02(1 GeV^) = 0.364 ± 0.126 respectively. These 
references favor a positive value for 02(1 GeV^) and the most recent one is done by Ref. 23| . 



-0.15 



0.17lo;i7 through a QCD light-cone sum rule analysis of the 



which shows that 02(1 GeV^' 

semi-leptonic i? — t- vr weak transition form factor based on the BABAR data on i? — t- tt/z/ 



28l | . The pion-photon transition form factor being involved only one pion DA maybe helpful 



to clarify the present situation. 

As argued in Ref. |l4| . the leading Fock state contributes to -^^^(0) only half and the 
remaining half should be come from the higher Fock states a.s Q^ -^ 0. And then both 
contributions from the leading Fock state and the higher Fock states are needed to get the 
correct 7r° — )■ 77 rate. In Ref. 29|], we have made such a comprehensive analysis of the 
pion-photon transition form factor in a smaller Q^ region, e.g. Q^ G [0, 10] GeV^, by taking 
both the valence quark and the non-valence quark contributions into consideration. It has 
been found that both the asymptotic-like and the CZ-like DAs can explain the CELLO 
and CLEO data [JJ, 12[ by setting proper parameters for the pion wavefunction. Then it 
shall be interesting to extend our previous analysis to higher Q^ region so as to determine 
which pionic behavior is more preferable for consistently explaining the CELLO, CLEO and 



^ Within the k± factorization approach, by keeping the transverse momentum dependence consistently and 
with the help of the Sudakov and threshold resummation, this end-point singularity may be cured to a 
certain degree, e.g. for pion-photon transition form factor 15| and for B — > light form factors 16[. 




FIG. 1: Typical diagrams that contribute to the pion-photon transition form factor F^^{Q'^), where 
x' = {1 — x). The rightmost shaded oval with a slant pattern stands for the strong interactions. 

BABAR data within the whole measured energy region. 

The paper is organized as follows. In Sec. II, we present the calculation technology to 
derive the valence and non-valence contributions to the pion-photon transition form factor. 
For such purpose, we construct a pion wavefunction model based on the BHL-prescription 
and present all the necessary formulae for discussion its properties. In Sec. Ill, we discuss 
what we can learn of the pionic leading Fock-state wavefunction/DA in comparison with 
CELLO, CLEO and BABAR experimental data. Some further discussion and comments 
are made in Sec. IV. 

II. CALCULATION TECHNOLOGY 

Generally, the pion-photon transition form factor 77* — > 7r° can be written as 

F.,{Q') = Fi^^\Q^) + Fi^^y\Q% (2) 

where F!;^\Q'^) is the usual valence quark part, FI^'^^^Q'^) stands for the non- valence quark 
part that is related to the higher Fock state of pion. The valence quark contribution F!;^\Q'^) 
dominates only as Q^ becomes very large. Fig.([T]) shows this point more clearly. F20{Q'^) 
comes from Fig.(IT^), which involves the direct annihilation of (gg)-pair into two photons, i.e. 
the leading Fock-state contribution that dominates the large Q^ contribution. F!;^^^\Q'^) 
comes from Fig.([T)D), in which one photon coupling 'inside' the pion wavefunction, i.e. strong 
interactions occur between the photon interactions that is related to the higher Fock states' 
contributions 30J. Under the light-cone perturbative QCD approach 4!], we can obtain 
the valence part FJi^^Q"^). While for the non- valence part F^i^^''(Q^), because of its non- 
perturbative nature, we shall construct a phenomenological model based on limiting behavior 



TABLE I: The explicit form of the spin-space wavefunction x {x,\i. 
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at Q^ — )■ and Q^ — )■ 00 to estimate it's contribution. 

Since the pion wavefunction is the key component of the pion-photon transition form 
factor, in the following subsections, we shall first make a discussion on its explicit form. 



A. Pion Wavefunction and Related DA 

Taking into account the Melosh rotation [3l| , the full form of the pion wavefunction can 

be written as 

vi/,,-(x,k^) = Y. x'^'^(a:,k^)vl/f_(x,k^), (3) 

A1A2 

where Ai and A2 are helicity states of the two constitute quarks, x'^^'^^(x, k_|_) stands for 
the spin-space wavefunction coming from the Wigner- Melosh rotation, ^'^^'^^(a;, k_|_) can be 



found in Refs. 



32 



35l |. whose explicit form is shown in TABjIl \l/^(x,kx) stands for the 



spatial wavefunction, and we adopt the factorized model to do our discussion, which is 
divided into a x-dependence part i^tt^x) and a k_L-dependence part. '{)tt{x) may or may not 
be the distribution amplitude, which depends on the explicit form of the k^-dependence 



part. Based on BHL prescription |l4j . 
written as 



32I-I 



the spatial wavefunction ^^^(x, k^) can be 



^f.(x,kx) = A(^^(a;)exp 



k^ +m 



2 n 



8/32x(l-x)J' ^^^ 

where the x-dependence part ip-n{x) can be expanded in Gegenbauer polynomials, and by 
keeping its first two terms, we obtain 



^J|(x, kx) = A (1 + 5 X Cl''^{2x - 1)) exp 



k^ +m^ 



(5) 



where the Gegenbauer polynomial C2 (2x 



■3/2/ 



8/32a;(l-x 
(3/2)[5(2x - If - 1]. The typical pa- 



rameter B determines the broadness of the wavefunction. The normalization constant A, 



the harmonic scale /3 and the hght constitute quark mass rriq are constrained by several 
reasonable constraints. The first is the conventional wavefunction normalization condition 

where /xq stands for some hadronic scale that is of order (9(1 GeV). The second is the 
constraint derived from vr*^ — )■ 77 decay amplitude [lj| 



0) ^ 



dx'^qq^X, kj 



^0 "- /tt 

Further more, rriq should be around the conventional adopted value 0.30 GeV. 
The leading Fock-state pion DA at the scale /iq takes the following form 



(7) 



0^(x,/io) 



2^/3 



t/2k| 



■'^nn(x,kl 



Substituting the wavefunction model ([3]), we obtain 

^/3AmP 



(8) 



Ax,f^l) = ^^^Zf^ \/^a^ (1 + ^ X cn^x - 1 



Erf 



\i 8/32x(l - x) 



Erf 



m^ 



\ 8/32x(l - x) 



(9) 



where the error function Erf (x) is defined as Erf (x) = -^ Jq e~* dt. Such a DA with 5 — )■ 
is asymptotic-like, and with the increment of B, it shall be broadened to a certain degree, 
e.g. when B ~ 0.6, it will be CZ-like with a fc_L-dependence factor which suppresses the 
end-point singularity. 

The pion DA at any scale Q^ can be derived from the initial DA 07r(x, /ig) through QCD 
evolution. The evolution equation up to order 0{as) takes the following form 



XX Q 



C; 



as{Q' 



[dy]V{x, y)(j)niy, f^o) - xx'(f)^{x, Q^ 



where [dy] = dydy'5{l — y — y'), (t)T,{x, Q"^) = 4>Tr{x, Q'^)/{xx') with x' = 1 — x, and 



(10) 



V{x,y) = 2CF 



xy'e{y - x) 6h,h2 + 



A 



y-x^ 



+ (1 f4 2) 



where (1 -H- 2) means that all the properties of the first constitute quark should be exchanged 
to that of second one, and vice versa. 5^^^^ = 1 when the two constitute quarks' helicities hi 
and /i2 are opposite and A0^(y, Q"^) = (pT^iy, Q^) — (pwix, Q^). With this evolution equation, 
we can take the evolution effects in calculating the pion-photon transition form factor. 



Moreover, a solution of Eq.(!TOl) in Gegenbauer expansion has been derived by Ref.|36|. 
which takes the following form 



0^(x, Q^) = 6xx' Y^ a„(/io) In 



Q' 



n=0 



A2 



cn2x-i] 



(11) 



where the Gegenbauer polynomials C^'^(2a; — 1) are eigenf unctions of V{x^y) and the cor- 
responding eigenvalues are the "non-singlet" anomalous dimensions 



In 



/3o 



' n+l 1 

1 + 4E 



25 



hiho 



k 



,^2'^ (n+l)(n + 2);' 

where /3o = 11 — 2n//3. The non-perturbative coefficients a„(/io) can be determined from 
the initial condition (j)T,{x,fil) by using the orthogonality relations for the Gegenbauer poly- 

i.e. 



nomials C^/^(2x 



«n(/io) 



J^dxMx,^l)C'J\2x-l) 



(12) 



J^ dxQx{l - x)[Cn^\2x - 1)]^' 

It should be noted that even though the model wavefunction (^ is constructed by using 
only the ffist two Gegenbauer terms in the longitudinal function (fTj{x), our present DA 
0^(x, /Xq) as shown by Eq.([n]) can be expanded in a full form of Gegenbauer series, i.e. both 
the leading and the higher Gegenbauer terms are there, whose corresponding Gegenbauer 
moments can be calculated with the help of Eq.( !T2|) . As will be shown in the following 
TABini the second Gegenbauer moment 02 (/ip) is close but not equal to the parameter B. 
This shows that the DA (p^^ is different from ^9^, which is due to the choice of the BHL- 
transverse momentum dependence and the consideration of all the helicity components' 
contributions. While by taking a simpler Gaussian-transverse momentum dependence and by 
taking only the usual helicity component into consideration, e.g. the transverse momentum 
dependence oc exp ( — r,^^(f_^;^ ] 10|, it leads to (J)t,{x,^1) = ^-wix), and a2(/^o) — -^• 



2ax{l-x) 



B. FiY^XQ"^) up to NLO 

Under the light-cone perturbative QCD approach ^, and by keeping the fc^-corrections 
in both the hard-scattering amplitude and the wavefunction, Ft^^{Q'^) has been calculated 



up to NLO 



14 



15 



29 



37 



39|. It is noted that for high helicity states (Ai + A2 = ±1), since 



their hard parts are proportional to the small current quark mass, we can safely neglect their 



contributions. As a combination of the LO part IJ, |29|, |37|, |38| and the NLO part 15|, |39| 



that keep the fc_L-dependence in the hard kernel, we can obtain the following formula after 
doing the integration over the azimuth angle, 



Fivm 



1 



1 p'Q' dx 



A-s/Stt'^ Jo Jo 



xQ2 



47r V 



/^/ 



xQ^ + kj 



+ 2 In X + 3 



vr 



(13) 



where [dx] = dxdx'S{l 



X 



x'), Cp = 4/3 and k^ = |k_|_|. fif stands for the factorization 



scale, and for convenience, we take fXf = Q 



dM- 



Here, without loss of generality, the usual 



assumption that the pion wavefunction depending on k^ through k\ only, i.e. "^qqlx, kj 
has been implicitly adopted. 



^qg{X,kl), 



F, 



(NV) 



7r7 



m 



As for FI^!^^\Q'^), due to its non-perturbative nature, it is hard to be calculated in any 



Q^ region. As stated in Ref. 



■•3 



around the region of Q^ ~ 0, we can treat the photon 



'inside' the pion wavefunction (nearly on-shell) as an external field that is approximately 
constant throughout the pion volume. And then, a fermion in a constant external field is 
modified only by a phase, i.e. Sa{x — y) = e~'"^'^y~^^'^Sp{x — y). Consequently, the lowest 
gg- wavefunction for the pion is modified only by a phase e~*^^'^, where y is the gg-separation. 
Transforming such phase into the momentum space and applying it to the wavefunction, we 
can obtain the two limiting behavior of F^!^'^\Q'^) at Q^ — )■ 0, which can be written as 

1 



i^ir(o)=Fi,^)(o) 



8V37r2 



dx'^qg{x,0±), 



(14) 



and 



d 



dQ 



Fir\Q' 



2 T7 



|Q2^0 



8^3: 



7r^ 



d 



1 r'Q' (^qq{x,kl] 



dXdrv I 



-A 



128v/3m27r2/32 

^2 1 



exp 



dQ^ Jo Jo V 2:2Q2 

C (l + Byi cTi2x - l)) -im' + Axx'(3^) ■ 
Jo ^ ' x' 

dx, (15) 



^(5'^xx' _ 

where x' = 1 — x. The above equation shows explicitly that at Q^ — )■ 0, the leading Fock 
state contributes to F^^(O) only half, i.e. F^^)(0) = F^^(0)/2. While by taking both the 



8 



valence and non- valence contributions into consideration, one can get the correct rate of the 



process n^ — )■ 77. 



Next, we construct a phenomenological model for FJ^!^'^\Q'^) by requiring it satisfy the 
above listed two limiting behavior at Q^ = and by assuming that it is power suppressed 
to F^{Q'^) in the limit Q^ —^ 00. For such purpose, we adopt the model constructed in 



Ref. 



29| 



Fi^^)(Q') (^^g2/^2)2' (16) 



where k = / — -g — ^my^, and a = iF^^(O). It is easy to find that FJ^^^\Q^) will 

be suppressed by 1/Q^ to F^^^(Q^) in the limit Q^ —^ 00. Then at large Q^ region, the 
non- valence Fock state part F^'^ \Q'^) shall give negligible contribution to the form factor. 
However it shall give sizable contribution at small Q^ region. 



D. Probability Pqq and Charged Mean Square Radius (r^+)'^^ 

After deriving the possible ranges for the parameters in the pion wavefunction, we shall 

meet the question that whether the resultant wavefunction and hence its DA is reasonable 

or not. In addition to the pion-photon transition form factor, the pion electromagnetic form 

actor F^+{Q^) also provides a platform for studying the properties of pion wavefunction 



34 



40 



43|. 



For such purpose, following the same procedure as described in detail in Ref. [3J], we derive 
a formula for the soft part contribution by taking all the helicity components' contribution 
to the pion electro-magnetic form factor. The general form for the soft part contribution 



can be written as 44 1 



/fj'j'c] \c I 
^^ E ^;-(^> k^> Ai)v^,,(x, k'^, A2), (17) 



dx(P\<i i_ 

where Q^ = qi and k'_L = k_L + (1 — 2;)q_L for the final state LC wavefunction when taking 
the Drell-Yan-West assignment. We can derive the probability for finding the lowest valence 
quark state Pqq and the charged mean square radius {r'^+Y^ from the limiting behavior of 

F^%(g2)atg2^o. 

Substituting the pion model wavefunction ([3]) and finishing the integration over k^ with 
the help of the Schwinger a— representation method, ^ = j^ /g°° a'^~^e~'^^da, Eq.( 1T7|) can 



be simplified as 



1 rca y^2 
dx / dX —, -TT7 

Jo 1287r2(l + A)3 



1 + Bx C2^\2x - 1) 



exp 



8m2(l + A)' + gV2(2 + A(4 + A)) 



2^/\2 



-Q'x'X 



32x'x/32(l + A) 
32x' X 13^ {I + \)- 



X 



'^ V32x/32(1 + A) 
QV^(2 + A(4 + A)) + Smlil + Xf] - h [^^^^^ Q'^'>^'] ' (18) 



where x' = 1 — x and J„ (n = 0, 1) stands for tlie modified Bessel function of tlie first 
kind. After taking tlie expansion in the small Q^ limit, we obtain the probability Pqq for the 
valence quark state, 

dxd'^k I 



<?<? 



KAQ' 



=0 



Ti^l'^^^-^^'^^^)!' 



1 ^oo y^2 

dx / dX — -— - 

io 167r2(l + A)2 



ml{l + X) + Axx' 13^ 



1 + Bx &J'^(2x- 1 



exp 



m^(l + A)' 



i^RI 



4:x'xf3 



(19) 



and the charged mean square radius (r^+)'^'^. 



2 \gg 



(r^) 



9Q2 



12 



2N^/ 



3^2(2 + 4A + A2)s 



(ix / dX- 

io 2567r2x/32(l + A)'^ 



1 + B X C'^^\2x - 1) 



X exp 



m^(l + A)' 



/^/32 



4:X'X(3 



ix'xl3'^ + ml{l + X) 



(20) 



In the above two equations, one may observe that the terms in the big parenthesis that are 
proportional to m^ come from the ordinal helicity components, while the remaining terms 
in the big parenthesis are from the higher helicity components. 



III. NUMERICAL RESULTS 



We adopt the NLO as{Q'^) to do the numerical calculation, i.e. 

47r 



as{Q') 



2^\nMQVKcDj[ 



(21) 



P,HQ^/KIcd) 

where /So = H — 2rT,//3 and /3i = 51 — 19r2//3. The value of nj varies with the energy scale 
and the value of Is.qcd is determined by requiring ^^(m^o) = 0.1184 [5|, i.e. Aqcd = 0.231 
GeV. 
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TABLE II: Pion wave function parameters under the condition of rriq = 0.30 GeV, and its proba- 
bility Pqq, charged mean radius w (r^+)'^^ (unit: fm) and the Gegenbauer moments a2,4,6(A'o)- 
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A{GeV-^) 


/3(GeV) 
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\/i^l+)'' 


«2(^o) 


04 (/io) 


aeilJ-o) 


0.00 


25.06 
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63.5% 


0.341 


0.027 


-0.027 


-0.016 
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21.71 


0.641 


60.0% 


0.358 


0.250 


-0.025 


-0.034 


0.30 


20.26 


0.668 


62.0% 


0.378 


0.362 


-0.018 


-0.041 


0.40 


18.91 


0.695 


66.1% 


0.401 


0.471 


-0.008 


-0.047 


0.60 


16.62 


0.745 


79.9% 


0.451 


0.679 


0.020 


-0.054 
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FIG. 2: Comparison of the pion DA model defined in Eq.([8]) with the asymptotic-form DA and 
the CZ-form DA, where B = 0.00, 0.30 and 0.60 respectively. 

A. Properties of Pion Wavefunction and Pion DA 



By taking /io = 1 GeV and mg = 0.30 GeV, we present the wavefunction parameters 
in Tabini which are determined by the mentioned constraints and by taking B = 0.00, 
0.20, 0.30, 0.40 and 0.60 respectively. The probabihty for the valence quark state Pqq, the 



charged mean radius J (r'^+Y^ (unit: fm) and the Gegenbauer moments a2,4,6(/^o) ^^^ ^^^ 
presented in TabJTTl For the case of i? = 0, a2,4,6(Aio) '^^'^ t)e safely neglected due to their 
smallness, then the corresponding DA is close to the asymptotic form as shown explicitly 
by Fig. ([2]). For a bigger B, it is found that a2(Aio) usually is quite larger than 04,6 (/^o); 



11 



which is consistent with our model wavefunction (^, where only the first two Gegenbauer 
terms are kept in (p{x). It is noted that by varying the parameter B within the region of 
~ [0.00, 0.60], the pion DA shall vary from aymptotic-like to CZ-like form. To show this point 
more clearly, we draw the pion DA defined in Eq.(|H]) in Fig. ([2]), where B = 0.00, 0.30 and 
0.60 respectively. As a comparison, we also present the conventional asymptotic-form DA, 
(f'Asix) = 6x(l — x) [J], and the CZ-form DA, (pczix) = 30x(l — x){2x — 1)^ ITJ. One may 



observe from Tab HTl that the value of {r'^+Y^ increases with the increment of 5, which runs 
within the region of [(0.341fm)^, (0.451fm)^] by varying B G [0.00,0.60]. These values are 
somewhat smaller than the measured pion charged radius {r'^)expt = (0.657=1=0.012 fm)^ 45 1 
and (0.641fm)^ [46|, but it is close to the value as suggested in Refs.|34J. I47ll48|. Such smaller 
(r^+)'^^ for the leading Fock-state wavefunction is reasonable, since the probability of leading 
Fock state Pqq is less than 1 and is about 60% — 80%. This confirms the necessity of taking 
the higher Fock-states into consideration to give full estimation of the pion electromagnetic 
form factor/pion-photon transition form factor, especially for lower Q^ regions. 

A naive pion wavefunction model has been suggested in Ref. 10|] to explain the new 
BABAR data J3[, which is constructed with a flat DA together with a Gaussian ansatz for 
the fc^-dependence, i.e. 



'^qq{x, kj 



47rV7r07r(a;) 



A;2 



VSxx'a """"^y 2axx')' ^^^^ 

where x' = 1 — x and 07r(x) = 1. With such a model (by set ting a = 0.53GeV^), it can be 
easily see that one can not derive the right behavior at Q^ — )■ [lO| , since following the similar 
steps as shown in Sec. II. D, it will lead to Pgg = /q ( g^^ J dx and (r^+)'^^ = /q [^t^] dx, 
both of which are divergent. This in some sense explains why such model wavefunction 
can explain the pion-photon transition form factor's large Q^ behavior (due to the large 
enhancement at the end-point region), but fails to explain the lower Q^ behavior. 

Next, it would be interesting to make a comparison with Brodsky and Teramond's holo- 



graphic model (BT model) with a quark mass effect |49l. l50l] for the pion DA. The BT model 



is predicted by using the anti-de Sitter / conformal field theory (AdS/CFT) correspondence 
and by using the soft-wall holographic model, whose explicit form is 50 1 

4>M{x,fil) 



C\J x{l — x) exp 
where k = 0.375 GeV |49|, m„ : 



1 / 2 2 

1 / mz m^ 



X I — X 



2k^ 
2 MeV and m^ = 5 MeV 



1 — exp 



2k'^x(1 ■ 



X] 



(23) 



50|. The factor C ~ 2.55, which 



can be determined by its normalization. It is found that for /xq ~ IGeV, the term involving 
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FIG. 3: Comparison of the pion DA model defined in Eq.Q with the Brodsky and Teramond's 
holographic model (BT model), where B = 0.00, 0.05, 0.10 and 0.15 respectively. 



/io gives quite small contribution and it can be safely neglected as is done by Ref.[50|- A 
comparison of our present pion DA model (IHl) with that of the BT model is presented in 
Fig.Q, where our model with B = 0.00, 0.05, 0.10 and 0.15 are presented by the circles, 
the dashed, the dash-dot and dotted lines respectively, and the BT model is drawn by a 
solid line. Since both models have similar transverse momentum behavior, it is natural to 
estimate that when setting B ^ 0.125, which corresponds to the same second Gegenbauer 
moment of BT model 02 (/^o) ~ 0.145, these two models shall lead to a similar behavior 
for the pion-photon transition form factor if calculating under the calculation technology as 
described in Sec. II. C ^. 



B. Pion-Photon Transition Form Factor 

First, we calculate the pion-photon transition form factor with the model wavefunction 
([3]) by taking rriq = 0.30 GeV and by varying B within the region of [0.00, 0.60]. The result 
is shown in Fig.(jlj), where the dash-dot line, the dotted line and the dashed line are for 



^ For such a calculation, one needs to be careful that the spin-space wavefunction for BT model should be 
changed accordingly, since to„ and nid are taken different values that is different to our present treatment 
of TO„ = md = niq. 
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FIG. 4: Q'^ Ftj^{Q'^) with the model wavefunction ^ by taking niq = 0.30 GeV and by varying 
B within the region of [0.00,0.60]. The dash-dot hne, the dotted hne and the dashed hne are for 
B = 0.00, B = 0.30 and B = 0.60 respectively. The solid line is the fitted curve ([1]) derived by 
BABAR [31]. 

B = 0.00, B = 0.30 and B = 0.60 respectively. As a comparison, we also present the 
BABAR fitted curve ([1]) in Fig.(jl]), which is shown by a solid line. For small energy region, 
Q^ ^ 15 GeV^, it is found that both asymptotic-like and CZ-like wavefunctions by adjusting 
the quark mass parameter can explain the CELLO, CLEO and BABAR experimental data, 
which agrees with the observation in Ref . 29|| . However, at large Q^ region, different behavior 
of DA (by varying B) shall lead to different limiting behavior. Typically, it is found that 
when Q^ — )■ oo, the Q'^F.n^{Q'^) for asymptotic-like wavefunction (with 5 = 0) tends to 
the usual limit 2/,^ ^ O.lSSGel^ [4]. So to explain the newly obtained BABAR data on 
high energy region, we need a broader DA other than the asymptotic one. It is found that 
with the increment of B (corresponding to a more broader DA as shown by Fig. ([2])), the 
estimated pion-transition form factor shall be more close to the BABAR data. Therefore 
the pion DA behavior will be determined If BABAR present measurement can be confirmed 
in the coming future. 

Second, we show how the leading valence quark and the non-valence quark contribute 
to the pion-photon transition form factor. We show the results for B = 0.60 in Fig. ([5]), 
where the solid line, the dotted line and the dashed line are for total contribution, the 
leading valence quark contribution and the non-valence quark contribution to the pion- 
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FIG. 5: Q^F^^{Q'^) with the model wavefunction ([3]) by taking m^ = 0.30 GeV and B = 0.60. The 
sohd Une, the dotted hne and the dashed Une are for total contribution, the leading valence quark 
contribution and the non-valence quark contribution to the form factor respectively. 

photon transition form factor respectively. Fig.(|5]) shows that the leading valence Fock- 
state contribution dominates the pion-photon transition from factor Q'^F^^^iQ'^) for large 
Q^ region, and the non-valence quark part is small in high Q^ region, but it shall provide 
sizable contribution to the low and intermediate energy regions. So one should consider the 
non- valence Fock states' contribution to Q'^Ft^^{{]) so as to explain the experimental data at 
both low and high Q^ region. 

Third, we make a discussion on the uncertainties caused by varying the value of uiq. For 
such purpose, we fix B to be 0.60 ^. From Fig. ([2]), one may observe that when B = 0.6, the 
DA is close to the CZ-form. As has been argued in Ref. 29|, for the case of CZ-like DA, in 
order to be consistent with the experimental data at low energy scale, mq should be within 
the region of 0.40lojo GeV. So we vary nig within the region of [0.30, 0.50] GeV to show the 
uncertainties. The results are shown in Fig. ([6]), where the solid line is for rriq = 0.40 GeV, 
and the shaded band shows its uncertainty. In the lower Q^ region, the upper edge of the 
band is for rriq = 0.50 GeV and the lower edge is for rriq = 0.30 GeV; while in the higher Q^ 
region, the upper edge of the band is for m^ = 0.30 GeV and the lower edge is for rriq = 0.50 
GeV. 



The case oi B = 0.30 is similar, only the range of niq should be shifted to 0.30_q j^q GeV [29j . 
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FIG. 6: Q'^Ft^^{Q'^) with the model wavefunction 1^ by fixing B = 0.60 and by varying niq within 
the region [0.30, 0.50] GeV. The soUd Une is for niq = 0.40 GeV, and the shaded band shows its 
uncertainty. 

IV. SUMMARY 



In the present paper, we have taken both the valence quark state's and the non-valence 
quark states' into consideration. The valence quark part is calculated up to NLO within 
the kT factorization approach and the non- valence quark part is estimated by a naive model 
based on its limiting behavior at both Q^ — > and Q"^ ^ oo. Our results show that (1) 
For Q^ ^ 15 GeV^, it is found that both asymptotic-like and more broader wavefunctions 
can explain the CELLO, CLEO and BABAR experimental data under reasonable choices 
of parameters. To be consistent with the new BABAR data at large Q^ region, we need a 
broader DA, i.e. the conventional adopted asymptotic DA should be broadened to a certain 
degree. (2) With suitable parameters for the pion model wavefunction that is constructed 
based on the BHL prescription, it is found that a more broader DA (with larger B) shall lead 
to a better agreement with the BABAR data. If BABAR confirms its present measurement, 
then pion DA should be broader, such as a CZ-like one with an improved behavior at the 
end point region. (3) The present adopted model of the pion wavefunction as shown by 



Eq.(j4l) shall present a basis for the application of the pQCD approach 



40 



43|. 
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